Measures of genetic diversification in somatic tissues at bulk and single-cell resolution.

Intra-tissue genetic heterogeneity is universal to both healthy and cancerous tissues. It emerges from the stochastic accumulation of somatic mutations throughout development and homeostasis. By combining population genetics theory and genomic information, genetic heterogeneity can be exploited to infer tissue organization and dynamics in vivo. However, many basic quantities, for example the dynamics of tissue-specific stem cells remain difficult to quantify precisely. Here, we show that single-cell and bulk sequencing data inform on different aspects of the underlying stochastic processes. Bulk-derived variant allele frequency spectra (VAF) show transitions from growing to constant stem cell populations with age in samples of healthy esophagus epithelium. Single-cell mutational burden distributions allow a sample size independent measure of mutation and proliferation rates. Mutation rates in adult hematopietic stem cells are higher compared to inferences during development, suggesting additional proliferation-independent effects. Furthermore, single-cell derived VAF spectra contain information on the number of tissue-specific stem cells. In hematopiesis, we find approximately 2 × 105 HSCs, if all stem cells divide symmetrically. However, the single-cell mutational burden distribution is over-dispersed compared to a model of Poisson distributed random mutations. A time-associated model of mutation accumulation with a constant rate alone cannot generate such a pattern. At least one additional source of stochasticity would be needed. Possible candidates for these processes may be occasional bursts of stem cell divisions, potentially in response to injury, or non-constant mutation rates either through environmental exposures or cell-intrinsic variation.

Mutation divergence over space in tumour expansion.

Mutation accumulation in tumour evolution is one major cause of intra-tumour heterogeneity (ITH), which often leads to drug resistance during treatment. Previous studies with multi-region sequencing have shown that mutation divergence among samples within the patient is common, and the importance of spatial sampling to obtain a complete picture in tumour measurements. However, quantitative comparisons of the relationship between mutation heterogeneity and tumour expansion modes, sampling distances as well as the sampling methods are still few. Here, we investigate how mutations diverge over space by varying the sampling distance and tumour expansion modes using individual-based simulations. We measure ITH by the Jaccard index between samples and quantify how ITH increases with sampling distance, the pattern of which holds in various sampling methods and sizes. We also compare the inferred mutation rates based on the distributions of variant allele frequencies under different tumour expansion modes and sampling sizes. In exponentially fast expanding tumours, a mutation rate can always be inferred for any sampling size. However, the accuracy compared with the true value decreases when the sampling size decreases, where small sampling sizes result in a high estimate of the mutation rate. In addition, such an inference becomes unreliable when the tumour expansion is slow, such as in surface growth.

The mutational landscape of the adult healthy parous and nulliparous human breast.

The accumulation of somatic mutations in healthy human tissues has been extensively characterized, but the mutational landscape of the healthy breast is still poorly understood. Our analysis of whole-genome sequencing shows that in line with other healthy organs, the healthy breast during the reproduction years accumulates mutations with age, with the rate of accumulation in the epithelium of 15.24 ± 5 mutations/year. Both epithelial and stromal compartments contain mutations in breast-specific driver genes, indicative of subsequent positive selection. Parity- and age-associated differences are evident in the mammary epithelium, partly explaining the observed difference in breast cancer risk amongst women of different childbearing age. Parity is associated with an age-dependent increase in the clone size of mutated epithelial cells, suggesting that older first-time mothers have a higher probability of accumulating oncogenic events in the epithelium compared to younger mothers or nulliparous women. In conclusion, we describe the reference genome of the healthy female human breast during reproductive years and provide evidence of how parity affects the genomic landscape of the mammary gland.

Coordinated inheritance of extrachromosomal DNA species in human cancer cells.

The chromosomal theory of inheritance has dominated human genetics, including cancer genetics. Genes on the same chromosome segregate together while genes on different chromosomes assort independently, providing a fundamental tenet of Mendelian inheritance. Extrachromosomal DNA (ecDNA) is a frequent event in cancer that drives oncogene amplification, dysregulated gene expression and intratumoral heterogeneity, including through random segregation during cell division. Distinct ecDNA sequences, herein termed ecDNA species, can co-exist to facilitate intermolecular cooperation in cancer cells. However, how multiple ecDNA species within a tumor cell are assorted and maintained across somatic cell generations to drive cancer cell evolution is not known. Here we show that cooperative ecDNA species can be coordinately inherited through mitotic co-segregation. Imaging and single-cell analyses show that multiple ecDNAs encoding distinct oncogenes co-occur and are correlated in copy number in human cancer cells. EcDNA species are coordinately segregated asymmetrically during mitosis, resulting in daughter cells with simultaneous copy number gains in multiple ecDNA species prior to any selection. Computational modeling reveals the quantitative principles of ecDNA co-segregation and co-selection, predicting their observed distributions in cancer cells. Finally, we show that coordinated inheritance of ecDNAs enables co-amplification of specialized ecDNAs containing only enhancer elements and guides therapeutic strategies to jointly deplete cooperating ecDNA oncogenes. Coordinated inheritance of ecDNAs confers stability to oncogene cooperation and novel gene regulatory circuits, allowing winning combinations of epigenetic states to be transmitted across cell generations.

Immune selection determines tumor antigenicity and influences response to checkpoint inhibitors.

In cancer, evolutionary forces select for clones that evade the immune system. Here we analyzed >10,000 primary tumors and 356 immune-checkpoint-treated metastases using immune dN/dS, the ratio of nonsynonymous to synonymous mutations in the immunopeptidome, to measure immune selection in cohorts and individuals. We classified tumors as immune edited when antigenic mutations were removed by negative selection and immune escaped when antigenicity was covered up by aberrant immune modulation. Only in immune-edited tumors was immune predation linked to CD8 T cell infiltration. Immune-escaped metastases experienced the best response to immunotherapy, whereas immune-edited patients did not benefit, suggesting a preexisting resistance mechanism. Similarly, in a longitudinal cohort, nivolumab treatment removes neoantigens exclusively in the immunopeptidome of nonimmune-edited patients, the group with the best overall survival response. Our work uses dN/dS to differentiate between immune-edited and immune-escaped tumors, measuring potential antigenicity and ultimately helping predict response to treatment.

Phenotypic plasticity and genetic control in colorectal cancer evolution.

Genetic and epigenetic variation, together with transcriptional plasticity, contribute to intratumour heterogeneity1. The interplay of these biological processes and their respective contributions to tumour evolution remain unknown. Here we show that intratumour genetic ancestry only infrequently affects gene expression traits and subclonal evolution in colorectal cancer (CRC). Using spatially resolved paired whole-genome and transcriptome sequencing, we find that the majority of intratumour variation in gene expression is not strongly heritable but rather 'plastic'. Somatic expression quantitative trait loci analysis identified a number of putative genetic controls of expression by cis-acting coding and non-coding mutations, the majority of which were clonal within a tumour, alongside frequent structural alterations. Consistently, computational inference on the spatial patterning of tumour phylogenies finds that a considerable proportion of CRCs did not show evidence of subclonal selection, with only a subset of putative genetic drivers associated with subclone expansions. Spatial intermixing of clones is common, with some tumours growing exponentially and others only at the periphery. Together, our data suggest that most genetic intratumour variation in CRC has no major phenotypic consequence and that transcriptional plasticity is, instead, widespread within a tumour.

The evolutionary dynamics of extrachromosomal DNA in human cancers.

Oncogene amplification on extrachromosomal DNA (ecDNA) is a common event, driving aggressive tumor growth, drug resistance and shorter survival. Currently, the impact of nonchromosomal oncogene inheritance-random identity by descent-is poorly understood. Also unclear is the impact of ecDNA on somatic variation and selection. Here integrating theoretical models of random segregation, unbiased image analysis, CRISPR-based ecDNA tagging with live-cell imaging and CRISPR-C, we demonstrate that random ecDNA inheritance results in extensive intratumoral ecDNA copy number heterogeneity and rapid adaptation to metabolic stress and targeted treatment. Observed ecDNAs benefit host cell survival or growth and can change within a single cell cycle. ecDNA inheritance can predict, a priori, some of the aggressive features of ecDNA-containing cancers. These properties are facilitated by the ability of ecDNA to rapidly adapt genomes in a way that is not possible through chromosomal oncogene amplification. These results show how the nonchromosomal random inheritance pattern of ecDNA contributes to poor outcomes for patients with cancer.

Quantification of spatial subclonal interactions enhancing the invasive phenotype of pediatric glioma.

Diffuse midline gliomas (DMGs) are highly aggressive, incurable childhood brain tumors. They present a clinical challenge due to many factors, including heterogeneity and diffuse infiltration, complicating disease management. Recent studies have described the existence of subclonal populations that may co-operate to drive pro-tumorigenic processes such as cellular invasion. However, a precise quantification of subclonal interactions is lacking, a problem that extends to other cancers. In this study, we combine spatial computational modeling of cellular interactions during invasion with co-evolution experiments of clonally disassembled patient-derived DMG cells. We design a Bayesian inference framework to quantify spatial subclonal interactions between molecular and phenotypically distinct lineages with different patterns of invasion. We show how this approach could discriminate genuine interactions, where one clone enhanced the invasive phenotype of another, from those apparently only due to the complex dynamics of spatially restricted growth. This study provides a framework for the quantification of subclonal interactions in DMG.

Evolution via somatic genetic variation in modular species.

Somatic genetic variation (SoGV) may play a consequential yet underappreciated role in long-lived, modular species among plants, animals, and fungi. Recent genomic data identified two levels of genetic heterogeneity, between cell lines and between modules, that are subject to multilevel selection. Because SoGV can transfer into gametes when germlines are sequestered late in ontogeny (plants, algae, and fungi and some basal animals), sexual and asexual processes provide interdependent routes of mutational input and impact the accumulation of genetic load and molecular evolution rates of the integrated asexual/sexual life cycle. Avenues for future research include possible fitness effects of SoGV, the identification and implications of multilevel selection, and modeling of asexual selective sweeps using approaches from tumor evolution.

Computational Approaches of Quasi-Static Compression Loading of SS316L Lattice Structures Made by Selective Laser Melting.

Additive manufacturing methods (AM) allow the production of complex-shaped lattice structures from a wide range of materials with enhanced mechanical properties, e.g., high strength to relative density ratio. These structures can be modified for various applications considering a transfer of a specific load or to absorb a precise amount of energy with the required deformation pattern. However, the structure design requires knowledge of the relationship between nonlinear material properties and lattice structure geometrical imperfections affected by manufacturing process parameters. A detailed analytical and numerical computational investigation must be done to better understand the behavior of lattice structures under mechanical loading. Different computational methods lead to different levels of result accuracy and reveal various deformational features. Therefore, this study focuses on a comparison of computational approaches using a quasi-static compression experiment of body-centered cubic (BCC) lattice structure manufactured of stainless steel 316L by selective laser melting technology. Models of geometry in numerical simulations are supplemented with geometrical imperfections that occur on the lattice structure's surface during the manufacturing process. They are related to the change of lattice struts cross-section size and actual shape. Results of the models supplemented with geometrical imperfections improved the accuracy of the calculations compared to the nominal geometry.

Evolutionary dynamics of neoantigens in growing tumors.

Cancers accumulate mutations that lead to neoantigens, novel peptides that elicit an immune response, and consequently undergo evolutionary selection. Here we establish how negative selection shapes the clonality of neoantigens in a growing cancer by constructing a mathematical model of neoantigen evolution. The model predicts that, without immune escape, tumor neoantigens are either clonal or at low frequency; hypermutated tumors can only establish after the evolution of immune escape. Moreover, the site frequency spectrum of somatic variants under negative selection appears more neutral as the strength of negative selection increases, which is consistent with classical neutral theory. These predictions are corroborated by the analysis of neoantigen frequencies and immune escape in exome and RNA sequencing data from 879 colon, stomach and endometrial cancers.

Subclonal reconstruction of tumors by using machine learning and population genetics.

Most cancer genomic data are generated from bulk samples composed of mixtures of cancer subpopulations, as well as normal cells. Subclonal reconstruction methods based on machine learning aim to separate those subpopulations in a sample and infer their evolutionary history. However, current approaches are entirely data driven and agnostic to evolutionary theory. We demonstrate that systematic errors occur in the analysis if evolution is not accounted for, and this is exacerbated with multi-sampling of the same tumor. We present a novel approach for model-based tumor subclonal reconstruction, called MOBSTER, which combines machine learning with theoretical population genetics. Using public whole-genome sequencing data from 2,606 samples from different cohorts, new data and synthetic validation, we show that this method is more robust and accurate than current techniques in single-sample, multiregion and longitudinal data. This approach minimizes the confounding factors of nonevolutionary methods, thus leading to more accurate recovery of the evolutionary history of human cancers.

Measuring the distribution of fitness effects in somatic evolution by combining clonal dynamics with dN/dS ratios.

The distribution of fitness effects (DFE) defines how new mutations spread through an evolving population. The ratio of non-synonymous to synonymous mutations (dN/dS) has become a popular method to detect selection in somatic cells. However the link, in somatic evolution, between dN/dS values and fitness coefficients is missing. Here we present a quantitative model of somatic evolutionary dynamics that determines the selective coefficients of individual driver mutations from dN/dS estimates. We then measure the DFE for somatic mutant clones in ostensibly normal oesophagus and skin. We reveal a broad distribution of fitness effects, with the largest fitness increases found for TP53 and NOTCH1 mutants (proliferative bias 1-5%). This study provides the theoretical link between dN/dS values and selective coefficients in somatic evolution, and measures the DFE of mutations in human tissues.

Exploiting evolutionary steering to induce collateral drug sensitivity in cancer.

Drug resistance mediated by clonal evolution is arguably the biggest problem in cancer therapy today. However, evolving resistance to one drug may come at a cost of decreased fecundity or increased sensitivity to another drug. These evolutionary trade-offs can be exploited using 'evolutionary steering' to control the tumour population and delay resistance. However, recapitulating cancer evolutionary dynamics experimentally remains challenging. Here, we present an approach for evolutionary steering based on a combination of single-cell barcoding, large populations of 108-109 cells grown without re-plating, longitudinal non-destructive monitoring of cancer clones, and mathematical modelling of tumour evolution. We demonstrate evolutionary steering in a lung cancer model, showing that it shifts the clonal composition of the tumour in our favour, leading to collateral sensitivity and proliferative costs. Genomic profiling revealed some of the mechanisms that drive evolved sensitivity. This approach allows modelling evolutionary steering strategies that can potentially control treatment resistance.

Measuring single cell divisions in human tissues from multi-region sequencing data.

Both normal tissue development and cancer growth are driven by a branching process of cell division and mutation accumulation that leads to intra-tissue genetic heterogeneity. However, quantifying somatic evolution in humans remains challenging. Here, we show that multi-sample genomic data from a single time point of normal and cancer tissues contains information on single-cell divisions. We present a new theoretical framework that, applied to whole-genome sequencing data of healthy tissue and cancer, allows inferring the mutation rate and the cell survival/death rate per division. On average, we found that cells accumulate 1.14 mutations per cell division in healthy haematopoiesis and 1.37 mutations per division in brain development. In both tissues, cell survival was maximal during early development. Analysis of 131 biopsies from 16 tumours showed 4 to 100 times increased mutation rates compared to healthy development and substantial inter-patient variation of cell survival/death rates.

The feedback between selection and demography shapes genomic diversity during coevolution.

Species interactions and coevolution are integral to ecological communities, but we lack empirical information on when and how these interactions generate and purge genetic diversity. Using genomic time series data from host-virus experiments, we found that coevolution occurs through consecutive selective sweeps in both species, with temporal consistency across replicates. Sweeps were accompanied by phenotypic change (resistance or infectivity increases) and expansions in population size. In the host, population expansion enabled rapid generation of genetic diversity in accordance with neutral processes. Viral molecular evolution was, in contrast, confined to few genes, all putative targets of selection. This study demonstrates that molecular evolution during species interactions is shaped by both eco-evolutionary feedback dynamics and interspecific differences in how genetic diversity is generated and maintained.

Spatially constrained tumour growth affects the patterns of clonal selection and neutral drift in cancer genomic data.

Quantification of the effect of spatial tumour sampling on the patterns of mutations detected in next-generation sequencing data is largely lacking. Here we use a spatial stochastic cellular automaton model of tumour growth that accounts for somatic mutations, selection, drift and spatial constraints, to simulate multi-region sequencing data derived from spatial sampling of a neoplasm. We show that the spatial structure of a solid cancer has a major impact on the detection of clonal selection and genetic drift from both bulk and single-cell sequencing data. Our results indicate that spatial constrains can introduce significant sampling biases when performing multi-region bulk sampling and that such bias becomes a major confounding factor for the measurement of the evolutionary dynamics of human tumours. We also propose a statistical inference framework that incorporates spatial effects within a growing tumour and so represents a further step forwards in the inference of evolutionary dynamics from genomic data. Our analysis shows that measuring cancer evolution using next-generation sequencing while accounting for the numerous confounding factors remains challenging. However, mechanistic model-based approaches have the potential to capture the sources of noise and better interpret the data.

How many samples are needed to infer truly clonal mutations from heterogenous tumours?

BACKGROUND: Modern cancer treatment strategies aim to target tumour specific genetic (or epigenetic) alterations. Treatment response improves if these alterations are clonal, i.e. present in all cancer cells within tumours. However, the identification of truly clonal alterations is impaired by the tremendous intra-tumour genetic heterogeneity and unavoidable sampling biases. METHODS: Here, we investigate the underlying causes of these spatial sampling biases and how the distribution and sizes of biopsies in sampling protocols can be optimised to minimize such biases. RESULTS: We find that in the ideal case, less than a handful of samples can be enough to infer truly clonal mutations. The frequency of the largest sub-clone at diagnosis is the main factor determining the accuracy of truncal mutation estimation in structured tumours. If the first sub-clone is dominating the tumour, higher spatial dispersion of samples and larger sample size can increase the accuracy of the estimation. In such an improved sampling scheme, fewer samples will enable the detection of truly clonal alterations with the same probability. CONCLUSIONS: Taking spatial tumour structure into account will decrease the probability to misclassify a sub-clonal mutation as clonal and promises better informed treatment decisions.

An explicit micro-FE approach to investigate the post-yield behaviour of trabecular bone under large deformations.

Homogenised finite element (FE) analyses are able to predict osteoporosis-related bone fractures and become useful for clinical applications. The predictions of FE analyses depend on the apparent, heterogeneous, anisotropic, elastic, and yield material properties, which are typically determined by implicit micro-FE (µFE) analyses of trabecular bone. The objective of this study is to explore an explicit µFE approach to determine the apparent post-yield behaviour of trabecular bone, beyond the elastic and yield properties. The material behaviour of bone tissue was described by elasto-plasticity with a von Mises yield criterion closed by a planar cap for positive hydrostatic stresses to distinguish the post-yield behaviour in tension and compression. Two ultimate strains for tension and compression were calibrated to trigger element deletion and reproduce damage of trabecular bone. A convergence analysis was undertaken to assess the role of the mesh. Thirteen load cases using periodicity-compatible mixed uniform boundary conditions were applied to three human trabecular bone samples of increasing volume fractions. The effect of densification in large strains was explored. The convergence study revealed a strong dependence of the apparent ultimate stresses and strains on element size. An apparent quadric strength surface for trabecular bone was successfully fitted in a normalised stress space. The effect of densification was reproduced and correlated well with former experimental results. This study demonstrates the potential of the explicit FE formulation and the element deletion technique to reproduce damage in trabecular bone using µFE analyses. The proper account of the mesh sensitivity remains challenging for practical computing times.